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Abstract 



lO ' For field equations of 4th order, following from a Lagrangian "Ricci 

o : 

scalar plus Weyl scalar" , it is shown (using methods of non-standard 
^ ■ analysis) that in a neighbourhood of Minkowski space there do not 

exist regular static spherically symmetric solutions. With that (be- 
sides the known local expansions about r = and r = oo resp.) for 
the first time a global statement on the existence of such solutions is 
^ ' given. Finally, this result will be discussed in connection with Ein- 

H : 

^ stein's particle programme. 

Fiir die Feldgleichungen 4. Ordnung, die aus einem Lagrangeaus- 
druck "Ricci-Skalar plus Weyl-Skalar" folgen, wird unter Zuhilfenahme 
von Methoden der Nicht-Standard- Analysis gezeigt, daB in einer Umge- 
bung des Minkowskiraumes keine statischen kugelsymmetrischen Losun- 
gen existieren. Damit wird erstmals neben den bekannten lokalen 
Entwicklungen um r = und r = oo eine globale Aussage iiber die 
Existenz solcher Losungen getroffen. Anschlief^end wird dies Ergebnis 
im Zusammenhang mit Einstein's Teilchenprogramm diskutiert. 
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1 Introduction 



General Relativity Theory starts from the Einstein- Hilbert Lagrangian Leh = 
R/2k with the Ricci scalar R which leads to Einstein's vacuum field equation 

— Q being of second order in the metrical tensor gi^. Its validity is proven 
with high accuracy in space-time regions where the curvature is small only. 
Therefore, the additional presence of a term being quadratic in the curvature 
is not excluded by the standard weak field experiments. 

Lagrangians with squared curvature have already been discussed by Weyl 
(1919), Bach (1921), and Einstein (1921). They were guided by ideas 
about conformal invariance, and Einstein (1921) proposed to look seriously 
to such alternatives. In Bach (1921) Lw = C/2k with the Weyl scalar 
2C = CijkiC'^^^^ was of special interest. Then it became quiet of them for 
a couple of decades because of the brilliant results of GRT but not at least 
because of mathematical difficulties. 

Recently, the interest in such equations has anewed with arguments com- 
ing from quantum gravity, cf. e.g. Treder (1975), BORZESZKOWSKI, 
Treder, Yourgrau (1978), Stelle (1978), Fiedler, Schimming (1983) 
and ref. cited there. In the large set of possible quadratic modifications a lin- 
ear combination L — L-er + I'^L^ enjoyed a special interest, cf. e.g. Treder 
(1977). The coupling constant I has to be a length for dimensional reasons 
and it has to be a small one to avoid confiicts with observations: One often 
takes Heisenberg length (= Compton wave length of the proton) = 1.3-10""^^ 
cm or Planck length 1.6 • 10~^^ cm. 

In the present paper wc consider this Lagrangian L in connection with 
Einstein's particle programme, see Einstein, Pauli (1943): One asks for 
spherically symmetric singularity- free asymptotically fiat solutions of the vac- 
uum field equations which shall be interpreted as particles, but this cannot 
be fulfilled within GRT itself; it is still a hope (cf. BORZESZKOWSKi (1981)) 
to realize it in such 4th order field equations. Two partial answers have al- 
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ready been given: Stelle (1978) showed the gravitational potential of the 
hnearized equations to be 

$(r) = — m/r + exp(— r//)c/r for r <^ / ; (1) 

he also gave an expansion series about r — 0, and Fiedler, Schimming 
(1983) proved its convergence and smoothness in a certain neighbourhood of 
r = 0, but said nothing about the convergence radius of this expansion. 

Now we want to join these two local expansions. To this end we write 
down the field equations in different but equivalent versions (set. 2), calculate 
some linearizations (set. 3) and prove the statement of the summary in sec- 
tion 4 which will be followed by a discussion of Einstein's particle programme 
in section 5. 



2 Notations and field equations 

We start from a space-time metric with signature (H ), the Riemann and 

Ricci tensor being defined by R^jki = ji,k ~ ■ ■ ■) ^ik = R'' ijk I'esp. The Weyl 
tensor C'^jki is the traceless part of the Riemann tensor and the Weyl scalar 
is given by 2C = CijkiC'^^^^. Light velocity is taken to be 1 and G — k/Stt is 
Newton's constant; g — det gij. 

Then we consider the Lagrangian 

C = y/^L = y/^{R + fC)/2K + £^at , (2) 

where Cmat is the matter Lagrangian. For writing down the corresponding 
field equation it is convenient to introduce the Bach tensor B^j, cf. Bach 
(1921) and WUNSCH (1976), beforehand. 

2 " V^Sg^^ - ^ -'"^ + 2^ ■ 

It holds Bj — 0, Bj.j — 0, Bij — Bji, and Bij conformally invariant of weight 
-1. 
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Variation of (2) with respect to Qij leads to the field equation 

Rij - \9ijR + l^Bij = KTij . (4) 

Now let us consider the vacuum case Tij = 0. The trace of eq. (4) then 
simply reads R = 0, and eq. (4) becomes equivalent to the simpler one 

Rij + l''Bij^O. (5) 

Writing k — and □• = g^^ one gets also the equivalent system 

R = 0; kRij + ORij = 2RiahjR'^^ + -gijRabR"'^ ■ (6) 

The static spherically symmetric line element in Schwarzschild coordinates 

reads 

ds^ = (1 + (3)e-^^dt^ - (1 + (3)-^dr^ - r^dQ^ , (7) 

where /3 and A depend on r only. The dot means differentiation with respect 
to r, and defining 

a = 2/3 - 2r/3 + 2rA(l + /?) + 2r^{X^ - A)(l + /3) + - 3A/3) , 

( = rd , Tji = ar~^ , rj2 = 3/3r~^ , 773 = (r~^ , 

the field equations (6) are just eqs. (2.15.a-c) and (2.16) of the paper of 
Fiedler, Schimming (1983). To avoid the products rrji we define a new 
independent variable x by r — le^ . Then rrji — rj^ — drji/dx and we obtain 

o^v'i + ^vi-m, (8) 

= (A; + rri3/6)r]2 + r^i + 773 + r^krji + r^riir]3/6 + rr]2r]3/2 + r^r]f/A , (9) 

= (1 + rr}2/SW3 - 2771 + 2A;772 + 2773 - r^krji - r^rjl/2 - r\rjs/6 + 2rrj2V3/^ ■ 

(10) 

Conversely, if the system eqs. (8-10) is solved, the metric can be obtained by 
/3 = r772/3, A = (2r772 + 2772 + rV)/(6 + 2r772), A(0) = 0. (11) 
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3 Linearizations 



Now we make some approximations to obtain the rough behaviour of the 
solutions. First, for r ^ oo, Stelle (1978) has shown the following: Lin- 
earization about Minkowski spaceQ leads (using our notations (1) and (7)) to 
A = and /? = 2$, where powers of $ are neglected and the term exp(r//)/r 
in $ must he suppressed because of asymptotical flatness. Further, a finite 
$(0) requires c = m in (1). Then $ is just the Bopp-Podolsky potential 
(stemming from 4th order electrodynamics), and it were Pechlaner, Sexl 
(1966) who proposed this form as representing a gravitational potential. 

Now we have $(r) = m (exp(— r/Z) — 1) /r = —m/l + rm/{2l^) + . . ., but 
also this finite potential gives rise to a singularity: the invariant RijR^^ ~ 
im? /r^ as r ^ 0, that means, the linearization makes no sense for this region. 
From this one can already expect that in a neighbourhood of Minkowski 
space no regular solutions exist. 

Second, for r ^ 0, Fiedler, Schimming (1983) proved that there exists 
a one-parameter family of solutions being singularity-free and analytical in a 
neighbourhood of r = 0. The parameter will be called e and can be defined 
as follows: neglecting the terms with r in eqs. (8-10) one obtains a linear 
system with constant coefficients possessing just a one-parameter family of 
solutions being regular at r = 0; it reads 

r/i = efc^r , r]2 = —5ekr , rj^ = Aek'^r . (12) 

(The factors fc" are chosen such that e becomes dimensionless.) Now one can 

take (12) as the first term of a power series rji = S„ a^^V" , and inserting 

this into eqs. (8-10) one iteratively obtains the coefficients aj"'\ For n even, 

vanishes: The r^-terms are k^{e/U + 10eV21), -k'^{e/2 - lOe^/3) and 

^That means, loosely speaking, the right hand side of eq. (6) will be neglected. Con- 
nected with this one may doubt the relevance of the term cr^^e^''/' in eq. (1) for m ^ 
and r — > oo, because it is small compared with the neglected terms. 
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k^{3e/7 + 20e^/7) resp. Furthermore, the r^"'~^-term is always a suitable 
power of k times a polynomial in e of the order < n. 

Up to the r^-terms the corresponding metric (7) reads 

ds^ = (1 + 5ker^/3)df - (1 + 5ker^/3)dr'^ - r^dil^ . (13) 

Third, we look for a linearization which holds uniformly for < r < oo. A 
glance at (11) and (7) shows that rji = gives Minkowski space, and therefore 
we neglect terms containing products of rji in (8-10). In other words, instead 
of the linearization used before we additionally retain the term r'^rji. Then 
we proceed as follows: from eqs. (8) and (10) we obtain 

^3 = + 3r7i and (14) 

_ -y'l - 5?7i + (e^^ - 4)r7i + l^e^'''r]i{i]i + i][/6) 

~ 2k + Ze^(r/; + 5r]i + Qr]i)/3 ' ^ ^ 

Inserting eqs. (14/15) into eq. (9) one obtains a third order equation for iji 

only0 whose linearization reads 

= < + 5r]'l + v[i2- e^") - Vii8 + 46^") . (16) 

The solution of eq. (16) which is bounded for x —oo reads| 

= 7 . [(3e-^^ + e-2^) sinh - Se"^^ cosh e^] . (17) 

A comparison with (12) gives 7 = 15el~^. 

^In MuLLER, Schmidt (1983) the same vacuum field equations were discussed for axi- 
ally symmetric Bianchi type I models. They possess a four-parameter group of isometrics, 
too, and the essential field equation is also a third order equation for one function. The 
difference is that in the present case spherical symmetry implies an explicit coordinate 
dependence. 

■^A second solution is 771 — — 12TOr^'* leading to the Schwarzschild solution (this solution 
solves both the full and the linearized equations in accordance with the fact that it makes 
zero both sides of eq. (6)), and the third one can be obtained from them up to quadrature 
by usual methods. 
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Now we insert this rji into eq. (15/11) and neglect again powers of rji, 
then the metric reads 
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P = be Ir ^ sinh(r/ ^) — cosh(r/ ^ 
{fz-'^ - 1) sinh(^r^) - Iz-'^ cosh(^r^)l dz . 



(18) 



This hnearization has (in the contrary to (1) and (12)) the following property: 
to each ro > and A > there exists an 5 > such that ioi —5 < e < 5 the 
relative error of the linearized solution (18) does not exceed A uniformly on 
the interval < r < Tq. 



4 Non-standard analysis 

In this section we will prove that in a neighbourhood of Minkowski space 
there do not exist any solutions. "Neighbourhood of Minkowski space" is 
in general a concept requiring additional explanations because of the large 
number of different topologies discussed in literature. To make the above 
statement sufficiently strong we apply a quite weak topology here: given 
a 5 > then all space-times being diffeomorphic to Minkowski space and 
fulfilling \RijE}^\ < Sk'^ form a neighbourhood about Minkowski space. 
In (13) we have at r = (independently of higher order terms in r) 

RijR'^ ^lQQk'^e'^/3 and Rqq ^ hke . (19) 

Therefore it holds: the one-parameter family of solutions being regular at 
r = is invariantly characterized by the real parameter e, and a necessary 
condition for it to be within a neighbourhood of Minkowski space is that e 
lies in a neighbourhood of 0. 

Now we suppose e to be an infinitesimal number, i.e., a positive number 
which is smaller than any positive real number. The mathematical theory 
deahng with such infinitesimals is called non-standard analysis, cf. Robin- 
son (1966). The clue is that one can handle non-standard numbers hke real 
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numbers. Further we need the so-called Permanence principle {— Robinson 
lemma): let ^(e) be an internal statement holding for all infinitesimals e. 
Then there exists a positive standard real 5 > such that A{e) holds for all 
e with < e < 5. The presumption "internal" says, roughly speaking, that 
in the formulation of A{.) the words "standard" and "infinitesimal" do not 
appear. For a more detailed explanation cf. the literature. 

This permanence principle shall be applied as follows: A[e) is the state- 
ment: "Take (12) as initial condition for eqs. (8-10) and calculate the corre- 
sponding metric (7/11). Then there exists an ro < Z such that RijR^^ > k'^ 
at r = Tq." 

Remark. At this point it is not essential whether tq is a standard or an 
(infinitely large) non-standard number. 

Proof oi A{e) for infinitesimals e with e 7^ 0: the difference of (18) in relation 
to the exact solution is of the order e^, i.e. the relative error is infinitesimally 
small. For increasing r, becomes arbitrary large, i.e. (1 + I3)~^ becomes 
small and RijR^^ increases to arbitrarily large values, cf. eq. (18). Now take 
To such that with metric (7/18) RijR^^ > 2k'^ holds. Then, for the exact 
solution, RijR^^ > k"^ holds at r = Tq, because their difference was shown to 
be infinitesimally small. 

Now the permanence principle tells us that there exists a positive standard 
real 6 > such that for all e with < |e| < 5 the corresponding exact solution 
has a point Tq such that at r = rg RijR'-' > k^ holds. But ''RijR^^ < fc^" is 
another necessary condition for a solution to lie within a neighbourhood of 
Minkowski space. 

Remark. Supposed this Tq is an infinitely large non-standard number, then by 
continuity arguments also a standard (finite) number with the same property 
exists. 

Now the statement is proved, but we have learned nothing about the ac- 
tual value of the number 6. Here, numerical calculations may help. They were 
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performed as follows: the power series for the functions rji were calculated up 
to the r^-term, then these functions taken at x — —A (i.e. r = 0.018 I) were 
used as initial conditions for a Runge-Kutta integration of eqs. (8-10). We 
got the following result: firstly, for e = ±10~^ and r < 101, the relative differ- 
ence between the linearized solution (17) and the numerical one is less than 
2 per thousand. Secondly, for < |e| < 1 the behaviour "/3 — >■ —oo ■ sgne" is 
confirmed. That means, the statement made above keeps valid at least for 

Remark. For large values e the power series for the rji converge very slowly, 
and therefore other methods would be necessary to decide about asymptotical 
flatness. 

5 Discussion — Einstein's particle programme 

Fourth order gravitational field equations could be taken as a field theoretical 
model of ordinary matter, the energy-momentum tensor of which is defined 

by 

kT*^ = - ^g.,R . (20) 

For our case one obtains at r = T*j to represent an ideal fiuid with the 
equation of state p* = /x*/3 and (cf. eq. (19)) 

lj,*{0) = bken-^ = 5e/87rG/2 . (21) 

Inserting |e| > 1 and Z < 1.3 • 10~^^ cm into eq. (21), we obtain 

|/i*(0)| > 1.5- 10^3 g cm-3. (22) 

Therefore it holds: if there exists a non-trivial static spherically symmetric 
asymptotically flat singularity-free solution of eq. (5) at all, then the corre- 
sponding particle would be a very massive one: its phenomenological energy 
density exceeds that of a neutron star by at least 40 orders of magnitude. 
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The resulting statement can be understood as follows: For a small curva- 
ture the 4th order corrections to Einstein's equations are small, too, and the 
situation should not be very different from that one we know from Einstein's 
theory. 

Now we want to refer to a problem concerning Schwarzschild coordinates: 
the transition from a general static spherically symmetric line element to 
Schwarzschild coordinates is possible only in the case that the function "in- 
variant surface of the sphere r = const, in dependence on its invariant radius" 
has not any stationary point. Here two standpoints are possible: either one 
takes this as a natural condition for a reasonable particle model or one allows 
coordinate singularities in (7) like "/3 > —1, and /? = — 1 at single points" 
which require a special care. (The discussion made above is not influenced 
by this.) 

The statement on the existence of solutions can be strengthened as fol- 
lows: Fiedler, Schimming (1983) proved that the solutions are analytical 
in a neighbourhood of r = 0. Further, the differential equation is an ana- 
lytical one and, therefore, in the subspace of singularity-free solutions they 
remain so in the limit r ^ oo. 

Then, there exists only a finite or countably infinite set of values e„ such 
that the corresponding solution becomes asymptotically fiat (the question, 
whether this set is empty or not, shall be subject of further investigation.); 
furthermore, the e„ have no finite accumulation point. That means, there 
exists at most a discrete spectrum of solutions. 

With respect to this fact, we remark the following: as one knows, Ein- 
stein's theory is a covariant one. But besides this symmetry, it is homotheti- 
cally invariant, too. That means, if ds^ is changed to e^^ds^ with constant Xi 
then the tensor Rij — \gijR remains unchanged, whereas the scalars R and 
''For /3 < — 1 one would obtain a cosmological model of Kantowski-Sachs type. Eqs. 
(2-11) remain unchanged for this case. 
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C will be divided by e^^ and e^'^ resp. From this it follows: with one solution 
of Einstein's vacuum equation one obtains by homothetical invariance just a 
one-parameter class of solutions. On the other hand, the sum R+l'^C has not 
such a symmetry and, therefore, one should not expect that a one-parameter 
family of solutions globally exists, and this is just in the scope of the particle 
programme where a definite particle's mass is wanted. 

Acknowledgement. Discussions with R. John, U. KASPERandR. SCHIM- 
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